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Abstract - We study the many-body physics in thin film topological band insulator, where the 
inter-edge Coulomb interaction can lead to an exciton condensation transition. We investigate 
the universality class of the exciton condensation quantum critical point. With different chemical 
potentials and interactions, the exciton condensation can belong to z = 2 mean field, or 3d XY, 
or Yukawa-Higgs universality classes. The interplay between exciton condensate and the time- 
reversal symmetry breaking is also discussed. Predictions of our work can be tested experimentally 
by tuning the chemical potentials on both surfaces of the thin film through gate voltage. We also 
show that all the analysis of the exciton condensate can be directly applied to a spin-triplet 
superconductor phase with attractive inter-edge interaction. 
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Quantum phases protected by topology have shown 
enormous interesting behaviors. Despite of unusual quan- 
tized responses to external fields [lip], topological phases 
usually manifest themselves with their stable edge states. 
For instance, the three dimensional topological band insu- 
lator (TBI) is characterized by its single Dirac cone edge 
state, with the Dirac point located at the time-reversal 
(T) invariant point in the edge Brillouin zone. These 
edge states were predicted theoretically, and also observed 
very successfully experimentally in materials such as BiSb, 
Bi 2 Se 3 and Bi 2 Te 3 [3j|7|[9]. It was understood that T is 
crucial to the stability of the edge states [7][8][lO][TT] , ba- 
sically because one cannot open up a T— invariant Dirac 
mass gap for a single edge. Enormous interests were de- 
voted to the T breaking at the edge states of TBI, includ- 
ing the cases with T broken by magnetic impurities and 
broken spontaneously due to interactions 
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In a thin film sample of TBI, since the two edges are 
close enough to interact with each other, T is no longer 
sufficient to protect the stability of the edge states i.e. it is 
allowed to open up a T invariant gap for both edge states. 
Recently, the gap at the edge states was indeed observed 
in experiments on thin film TBI, when the thickness of the 
film is small 
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It was proposed that even without 
direct inter-edge tunnelling, the local Coulomb interaction 
can also gap out the edge states through exciton conden- 
se, a particle-hole pair bound state across the 



specific biased gate voltage is applied to two edges, where 
there is a "nesting" between two Fermi surfaces, and the 
"exciton" susceptibility diverges [24]. Lately it was pro- 
posed that large dielectric constants of these materials in- 
creases the layer separation range over which the inter 
surface coherence survives 21 . In our current paper we 



hope to go beyond the mean-field consideration in Ref. 24 
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sation 

thin film condenses. This effect is most prominent when a 



and study the critical properties of the exciton condensate 
physics. 



We illustrate our main results of this paper in the phase 
diagram (Fig. [T]) plotted against the chemical potentials 
of the two edges, which we assume can be separately 
tuned with gate voltages on both edges. The color in this 
phase diagram denotes the critical Coulomb interaction 
U c required to drive the exciton condensate, and at the 
line ^1+^2 = except for the origin, U c = due to 
the divergence of exciton susceptibility. In most area of 
this phase diagram, the quantum critical point (QCP) be- 
longs to the z = 2 mean field universality class. At the 
origin /ii = /i2 = 0, the transition is described by the 
Yukawa-Higgs theory; and at the special line /ii = A*2 7^ 0, 
the transition belongs to the 3d XY universality class. 
The same phase diagram and universality classes apply to 
the superconducting transition with attractive inter-edge 
"Coulomb" interaction. 
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Fig. 1: The phase diagram of Eq. [T] The dark/bright colors 
denote the low/high critical Coulomb interaction. 



The effective Hamiltonian of the surface states reads 



H= ^2 ipj(via-jp- fJ-i)ipi + Un 1 n 2 , 



(1) 
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where v\ = (— l) l+1 Vf because the two edges of the TBI 
have opposite helicities. This Hamiltonian clearly has 
timc-rcvcrsal symmetry T : tpi — > ia y ipi, k — > —k. There 
is also a inversion symmetry X : ipi •<-> ife, k — > — k when 
Hi = H2, and the inversion is also a generic symmetry 
of materials such as Bi 2 Te3 and Bi 2 Se3. Gate voltages 
determine the relative chemical potentials, The sec- 
ond term of Eq. [I] describes the short ranged inter-edge 
Coulomb interaction. We assume the screening length of 
Coulomb interaction is always larger than the thickness 
of the think film, therefore the inter-edge Coulomb inter- 
action is still important even when the direct inter-edge 
electron tunnelling is ignorable. The intra-edge Coulomb 
interaction is also tentatively ignored in this Hamiltonian, 
its effects will be discussed later. 

Without inter-edge tunnelling, electron number in 
each layer is independently conserved. This enlarged 
U(l) x U(l) symmetry is spontaneously broken down to 
the diagonal U(l) symmetry by exciton condensation with 
order parameter <f> ~ U{ijj\tl>2), which also lowers the en- 
ergy of the system. Then a mean field Hamiltonian can 
be obtained from the Hubbard-Stratonovich transforma- 
tion of the original Hamiltonian Eq. [T| 24 : 



H 
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The complex order parameter <j> is invariant under T, but 
becomes its complex conjugate under inversion I. Clearly 
there are also other mean field channels of the Coulomb in- 
teraction such as {il)\a z 4 : 2), but the exciton order U(jp\if>2) 
has the lowest mean field energy because it opens up a 
Dirac mass gap, therefore we will focus on this order pa- 
rameter in our current work. 

After integrating out the fermions, the static and uni- 
form renormalization to the mean field Hamiltonian of d> 



Fig. 2: a, one-loop Feynman diagram for renormalized Boson 
Lagrangian, b, Feynman diagram for the fermion self-energy. 
The dashed lines can denote either the exciton (j> or exciton 
bilinear \4>\ 2 . The solid line denotes the fermion propagator. 



is C e ff = {jj — Xo)|0| 2 , Xo is the susceptibility of order pa- 
rameter <fi which can be evaluated from Boson self-energy 
loop diagram Fig. [2j 

Xo =5^(0,0)=- f Tr(Gi(k, iuj)G 2 (k, iu)), (2) 

J iu),k 

where Gi is Green's function of fermion. This integral 
increases linearly with the ultraviolet cut-off A i.e. it is 
not just a Fermi surface effect. For instance, when fix = 
fJ-2 — Mi Xo ~ A — \fi\. Had we included the other mean 
field order parameter ij}\a z il)2, since the susceptibility of 
this order parameter only comes from the Fermi surface, it 
would never beat the order parameter <f> ~ U(ip\ip2) under 
consideration, as long as the size of the Fermi surface is 
small compared with the UV cut-off. 

If we fix yit2, the critical Coulomb interaction U c is plot- 
ted in Fig. [3] As we mentioned, U c itself is cut-off de- 
pendent. However, the relative difference between U c is 
UV cut-off independent, hence we can still compare U c 
with no ambiguity. As we can see, at \i\ = — /i 2 ^ 0, 
the critical Coulomb interaction is zero, due to the fact 
that the exciton susceptibility diverges logarithmically at 
this point. The logarithmic divergence is a consequence of 
the Fermi surface nesting between the two edges [24] . At 
A*i = A*2 7^ 0, although the two Fermi surfaces still have 
the same size, the wavefunction overlap (V'ljjIV^fe) = 

for any momentum fc at the Fermi surface, therefore the 
susceptibility is not divergent. This matrix element sup- 
pression is due to the opposite helicity of the two edges, 
and as we will see, this suppression will also affect the 
dynamics and universality class of QCP. 

Now let us go beyond the mean field formalism and 
move on to the universality class of the QCPs. Without 
fermions, the exciton condensation would certainly belong 
to the 3d XY universality class, while coupling to Fermions 
will likely modify the universality class. Let us first discuss 
the case with fix = /i 2 = 0. At this point, after redefin- 
ing tp2 0z?/> 2 , the phase transition is described by the 
following Lagrangian: 

+ |^| 2 +^|V^| 2 +r|^| 2 + ^| 4 + --- 
rf=(rf, rf>), $= (Re[0], In#]). 
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Fig. 3: The plot of critical Coulomb interaction U c against 
chemical potential Hl> with fixed fi2 > 0. Critical Coulomb 
interaction is measured in units of the U c with = — 2. 

U c is zero at [i\ = —fi2 due to the logarithmically divergent 
exciton susceptibility. Everywhere but fii = ±/Lt2, the QCP 
have z — 2 dynamical exponents, with mean field universality 
class. The special point fj,l = /i2 7^ (solid circle) has 2 = 1 
dynamical scaling, and it belongs to the 3d XY universality 
class. 



Here 7^ = (a z , <r x , <t v ), rf and rj y are two Pauli matri- 
ces that mix the two edges. This model becomes pre- 
cisely the Higgs- Yukawa model which describes the Chi- 
ral symmetry breaking of Dirac fermion, at least when 
Vf — va. The transition of <p is not 3d XY transition be- 
cause the coupling A is relevant at the 3d XY fixed point, 
based on the well-known scaling dimensions [tp] = 1, and 
[$\ = (d - 2)/2 + 77/2 = 0.519 at the 3d XY fixed point 
1 22 . The critical exponents of this transition with large 
N have been calculated by means of 1/N and e = 4 — d 
and a second order transition with 
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expansions 

non- Wilson-Fisher universality class was found. Therefore 
we conclude that the transition is still second order, with 
different universality class from the 3d XY transition. 

If fii H2, due to the mismatch of the size of the Fermi 
surface, the loop diagram Fig. [2] will not lead to singu- 
lar behavior for Boson self-energy at low frequency and 
small momentum. However, the fermi surface mismatch 
breaks the inversion symmetry of the Hamiltonian Eq. [T] 
therefore the following term is allowed in the Lagrangian: 

d = h(i<fd T <i> y - i(j) v d T <t> x ) ~ h(/)*d T (f>, fin/a-m. (3) 

C\ leads to a z — 2 dynamical exponent, which is analo- 
gous to the Mott insulator (MI) -superfluid transition in 
Bosc Hubbard model away from the tip of MI lobe [28] , 
and also the XY magnetic transition in magnetic field. In 
two spatial dimension, the z — 2 transition is a mean field 
transition with marginally irrelevant perturbations. 

If fix = fi 2 7^ 0) then the exciton condensation is sim- 
ilar to a ferromagnetic transition in Fermi liquid, which 
usually has z = 3 over damped quantum critical mode 
29 30 . However, in our case there is the matrix element 



suppression effect mentioned before, namely: 

\M k , k+q \ 2 = |(/c + (Z ,l|^ +g ^ fe |fc,2)| 2 
= S m 2 (0 k+q -9 k )~q 2 /k 2 f . 



(4) 



This matrix element suppression strongly affects the low 
energy dynamics of the quantum critical fluctuations. For 
instance, the damping rate of the quantum critical modes 
due to particle-hole excitations can be calculated through 
the Boson self-energy diagram Fig. [2h: 



Im[E (w,g)] - g 2 



d 2 k 

(2^)2 I/( e *+*) 



/(<*)] 



x5(\uj\ - e k+q + e k )\M k . k+q \ 2 - g 



2 Mg 

Vfk 2 / 



(5) 



This term will not lead to over-damped z — 3 quantum 



critical modes. The same effect was noticed in Ref. 15 
in the context of T breaking at the edge state of TBI, 
and following the argument of Ref. 15 , we can conclude 



that this transition still belongs to the 3d XY universality 
class even though the order parameter <f> couples linearly 
to Fermi surface in Eq. [2] 

The quantum critical modes also modify the Fermion's 
self-energy. Using the Feynman diagram Fig. [2f), the 
Fermion self-energy reads 



E(zw) 



d 2 kdeG(ie + iw, k)(</>. s </>. z)\M , k \ 



Here we use the full correlation at the 3d XY fixed point: 
(<?i e k^-ie -%) ~ ( g2 + v 1^ 2 ) ~ 1+r1 ^ 2 , rj is the anomalous 
dimension of the order parameter <j> at the 3d XY fixed 
point. The leading order contribution to the imaginary 
part of Fermion self-energy reads 



£(w)" - (A) 2 M 2 +"sgn[ W ] « M 



(7) 



Therefore the Fermion self-energy correction is always 
dominated by the linear frequency term of the free elec- 
tron propagator, the linear coupling A does not destroy 
the Landau quasiparticles. 

In addition to the linear coupling Eq. [2] another 
quadratic interaction is also allowed by symmetry, but was 
omitted in the mean field Hamiltonian Eq. [2j 



£'~a'5>;<m<?) 2 



(8) 



When Hi = A*2 = 0, this term is clearly irrelevant based 
on straightforward power-counting. With finite Fermi sur- 
faces, after integrating out the fermions, this quadratic 
interaction will induce the following four-body interaction 
between d>: 



£2 = u(0)L,(f- 



(9) 
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Note that the term with |w| comes from the Fermi surface 
effects. We want to estimate the scaling dimension of u 
at the 3d XY fixed point. Since the scaling dimension 
[{<j>) 2 \ = 3— 1/V, the scaling dimension of u is [u] = 2/v— 3, 
here v is the standard critical exponent defined as £ ~ r~ v . 
Therefore as long as v > 2/3, u is irrelevant. This criterion 
is indeed satisfied according to the well-known exponents 
of 3d O(N) universality class 22 . £ 2 also exists at the 
z — 2 QCP with fix 7^ fi2 discussed before, however, since 
there [lj] = 2[q] = 2, this term has a high scaling dimension 
and is clearly irrelevant. A similar analysis about the \uj\/q 
term was first made in Ref. [3l] . 

Again we can evaluate effect of the coupling C on the 
Fermion self-energy. The leading order correction can 
again be calculated through diagram Fig. |2j>, while now 
the dashed lines are correlation functions (6 2 -6 1 -) ~ 

ie,k —ie, — k 

(e 2 + v^k 2 )2~". With finite Fermi surface, the leading 
order contribution to the imaginary part of fermion self- 
energy reads 



E(w)" ~ (A') 2 M 5 - 2/l/ sgnM « 



(10) 



Hence this quadratic coupling C does not destroy the Lan- 
dau quasiparticle at the QCP either. 

Now let us turn on an extra intra-edge Coulomb inter- 
action in Eq. [T] 



(11) 



Thi s term will favor to develop magnetization on each edge 
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Since Vm A n h V(ipja z ipi) 2 /2, the Hubbard- 

Stratonovich transformation can give us mean field or- 
der parameter <£>i ~ (ip^a z ip) and <J>2 ~ (tp^ a z r] z ip) with 
Ising symmetry. Without exciton condensate in the back- 
ground, these two Ising order parameters $1 and $2 are 
degenerate at the mean field level. $1 breaks only T, while 
$ 2 breaks both T and X. In the background of exciton 
condensate, $2 has lower fermion mean field energy be- 
cause the exciton order parameters anticommute with $2, 
hence the exciton condensate favors to have a I breaking 
magnetization. 

In Ref. [24], it was shown that when [i\ + (1% = 0, at 
the vortex core of the exciton condensate order parameter 
there is a Fermion zero mode, which carries charge \ . This 
zero mode is protected by the symmetry of Hamiltonian 
Eq. [1] 72 tf*72 = H, and 72 = ia v rf . This symmetry 
guarantees that the spectrum is symmetric with E = 0, 
and it is valid even with the presence of exciton vortex. 
With nonzero $2, this symmetry is broken, and there is 
no longer a zero mode at vortex core. By contrast, if the 
system develops magnetization $1, there is still a vortex 
core Fermion mode at precisely zero energy. 

Just like the exciton order parameter <j>, the order pa- 
rameter $ Q also couples to the Fermions both linearly and 
quadratically. Due to the same matrix element suppres- 
sion effect as in Eq. [5j the linear coupling does not lead 
to singular corrections to the QCP of $ a . However, using 



the similar argument as that below Eq. [9j the quadratic 
coupling C ~ \J2i Wi^ii^a) 2 will lead to a relevant per- 
turbation at the 3d Ising fixed point as long as one of the 
edges has a finite Fermi surface, due to the fact v < 2/3 



at the 3d Ising universality class 22 



In addition to the Fermi surface, the Goldstone mode of 
the exciton condensate couples to order parameter $ a as 
well, if the QCP of $ a occurs in a background of exciton 
condensate. In the case without X ^ /12), the lowest 
order coupling reads C" ~ A"(<9 T 6>)($ a ) 2 . 8 is the phase 
angle of the exciton condensate: <p ~ e l6 '. With C", after 
integrating out the Goldstone mode 9, a singular term is 
induced for 3> a : 



£3 = «3($a)L,g 



2 + vlq 2 



(*a) 2 _ iw , 



(12) 



To determine the scaling dimension of this term, we again 
have to compare v of 3d Ising transition and 2/3: since 
v < 2/3, this coupling £3 is also relevant at the 3d Ising 
universality class. This relevant perturbation exists even 
when the Fermions are fully gapped out by the exciton 
condensate. In the case with I, the coupling C" is for- 
bidden, since 9^—6 under X, In this case the coupling 
between $ a and the exciton Goldstone mode will occur at 
higher order, hence no relevant perturbation is induced at 
the 3d Ising universality class. 

Although the exciton is charge neutral, its transport ef- 



fect has been verified in bilayer quantum Hall system 35 



by measuring the tunnelling conductance between the two 
layers 34 . A similar measurement can in principle be 
carried out in the thin film topological insulator. The 
exciton condensate will lead to a sharp peak of the inter- 
surface tunnelling conductance. Inside the exciton con- 
densate phase, the condensate will be destroyed by the 
thermal fluctuation through a Kosterlitz-Thouless transi- 
tion at finite temperature. The scaling between the critical 
temperature T c of this KT transition and the tuning pa- 
rameter r depends on the universality class of the QCP, 
and r can be taken as the interaction U — U c . For exam- 
ple, for the z = 2 mean field transition, T c ~ \r\; while 
for the 3D XY transition in the phase diagram (Fig. [T]), 
T c ~ \r\ zv ~ |r| 2 / 3 . Thus different quantum critical be- 
haviors can be measured through T c . The interaction U 
between the two surfaces can be tuned by changing the 
thickness of the thin film sample. 

If U in Eq. [T] is attractive instead of repulsive, then 
the system favors to have superconductor pairing. Af- 
ter a particle-hole transformation for ip2 v x ip2' 
both U and /x 2 change sign, while all the other terms 
of the Hamiltonian remain unchanged. The most ener- 
getically favored pairing state is ipi^^, because after 
particle-hole transformation this pairing becomes the ex- 
citon condensate (j>. Therefore all the analysis of the QCP 
and Goldstone mode about this superconducting state can 
be obtained by particle-hole transformation of the exciton 
case. For instance, at chemical potential \i\ = fj>2, there 
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is a logarithmic divergence of pairing susceptibility, while 
at /xi = —fi2 there is a matrix element suppression at 
the interaction vertex between Cooper pair and fermions. 
The pairing ip\a x ^2 is a spin triplet pairing with total 
S* = 0, and the vortex core of this superconductor carries 
a fermion zero mode when fi± = 

In summary, we have studied the exciton condensation 
phase transition and its quantum critical properties in a 
phase diagram with edge dependent chemical potentials 
and Coulomb interaction. Interplay between exciton con- 
densate and other order parameters are also discussed. In 
addition to the TBI materials that are currently under in- 
tensive experimental studies, we expect our formalism to 
be applicable to TBI with strong correlation, for instance 
the materials with 5d electrons 
proposed recently. 
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